ON BAXTER'S DIFFERENCE SYSTEMS 



J. S. GERONIMO, L. GOLINSKII 

Abstract. We study the asymptotics of solutions of a difference system 
introduced by Baxter by using the general method for the asymptotic 
representation of such solutions due to Benzaid and Lutz. Some results 
of Tauberian type are obtained in the case when the spectral parameter 
belongs to the unit circle. 



1. Introduction 

In this paper we study the asymptotics of certain solutions to a system 
of difference equations introducded by Baxter [1, 2, 3]. Baxter himself was 
interested in generalizing some results that had been obtained by Szego, 
Geronimus, and Verblunsky on the difference system satisfied by polynomi- 
als orthogonal on the unit circle. This system has been of much interest 
since it appears in the study by Miller et al [13] and Gesztesy et al [9, 10] 
of solutions of the Ablowitz-Ladik equations and by Geronimo et.al [6] who 
generalized the results of Geronimo and Johnson [7] (see also [15, 16]). 



2. The Baxter equations 

Let [i be a complex Borel measure on the unit circle T of bounded to- 
tal variation. We assume throughout that the n-th Toeplitz determinant 
associated with fj, 



D n (//) = det||/i i _ i ||?Ti o ^0, 
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for all n = 1, 2, Construct polynomials <j) n and ip n of degree n in z and 

K respectively, by 
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such that, 

/ ^„(C)C~ fc ^= / ^ n (C)C fe dM = Ah-i(m)<U, A; = 0,1,...,: 
jt jt 

Pick £ n from £^ = (D n D n+ i)~ l [l n is determined up to a ± sign) and put 
0n (z) := ^n4>n(z) = K„2 n + . . . , 

Y^O) := Znlpniz) = H n Z~ n + ..., K n = £ n D n (fJ,) ^ 0. 

So 

/ 4(C)C^M = / MCK k d^ = — , fc = 0,l,...,n, 

JT JT K n 

and {(/'njV'n} form a bi-orthogonal sequence 

/ 0n(C)^m(C)^M = <*ra,m- 
JT 

The polynomials will be unique once a determination is made for their lead- 
ing coefficients. 

The Baxter parameters are introduced in the following way 



It is easy to see that 

Pn = -tl^-l(^) 
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Similarly, 
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Next, we have 
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On the other hand, the Silvester identity applied to Z) n +i(/z) gives 
Ah-i(aOAi-i(aO = D 2 n {n) - £> n (C/i)A»(CV), 

so 



1 - a n (3 n 



Since < = t 2 n D 2 n {n) = D n (ji)/D n+1 (ji), then 



1 - a n /3 n 



v n-l 



*n 2 = I[(l-aiPj)- 

3=1 



(2.1) 



We define two sequences {u n } and {v n } of polynomials of degree n in z 
and i respectively, by 



u n (z) := K n z n i) n (z) = 
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(2.2) 



(2.3) 



It is easy to check that f3 n u n and n n (p n — n n -iz(f) n -i, both polynomials in 
z of degree < n, have the same Fourier coefficients with respect to fi. The 
same goes for <y. n v n and K n ifj n - K n -\z 1 ip n -i. Hence 

P n U n {z) = K n (j) n (z) - K n -lZ<f) n -l(z), 
a n V n {z) = Kn^n(z) ~ K«-l^Vn-l(4 

In the matrix form 

V 1 -z~^ n 
-a n 1 



$n-l(z) 



Kn—l 



^n(z), 



4>n{z) 




<t>n{z) 


Z n ^ n {z)_ 




K- 1 U n {z)_ 



Taking inverse matrices and shifting indices we have by (2.1) 



$n+l(*) 



Kn+1 



Z fin+1 

a n+ iz 1 



$n(z). 

*„(z), = 



Likewise we obtain 

Inspection of the above equations shows that 

tf n (z;a,0) = <S> n {z-\p,a) 



1pn(z) 
K^Vniz) 



(2.4) 

(2.5) 
(2.6) 



(since </> = Y>o = !)■ 

In the case when j5 n = a n , \a n \ < 1 for all n, the equations of Baxter 
reduce to those satisfied by polynomial orthogonal with respect to some 
positive measure supported on the unit circle. 

3. The Benzaid and Lutz method 

As mentioned in the introduction, Benzaid and Lutz [4] developed meth- 
ods for the asymptotic representation for solutions of difference equations. 
We will be interested in the case when the coefficients in equation (2.4) 
satisfy 

oo oo 

^|a n | 2 <oo, ^|/5 n | 2 <oo, and a n (3 n / 1 Vn > 1. (3.1) 



n=l 



n=l 



If in equation (2.4) we switch to the monic polynomials <I> n = K n 1( I> n , it 
becomes 



®n+l{z) = 

U{z,n) = 



a n+ iz 1 

Pn+l 
Oin+lZ 



<D n (z) = [A(z) + U(z,n)\ $ n (z), 
, A(z) 



z 
1 



(3.2) 



Q(z,n) 



The unperturbed system is Q^+ii 2 ) = M z )^n( z )i an d following Benzaid 
and Lutz we choose a Q(z,n) such that diag(Q(z, n)) = satisfying the 
equation 

U(n) + A{z)Q{z, n) - Q(z, n + l)A(z) = 0. (3.3) 
It is easy to check that Q in (3.3) can be taken as 

qn(z,n) 
_q 2 i(z,n) 

with 



(3.4) 



We could have defined a new variable X from & n (z) = (I + Q(z,n))X(z,n), 
and so 

X(z, n + 1) = (I + Q(z, n + l))" 1 (A(z) + U(z, n)){I + Q(z, n))X(z, n). 

However in order not to introduce a denominator when computing inverses 
(/ + Q)^ 1 we will carry out the Benzaid and Lutz procedure in two steps 
using nilpotent matrices 



qi2{z,n) 
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= Y J p k z n ~ k = 
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-j 4. , 
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oo 

- E akzk 
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-n __ 


OO 

^a n+j z 3 

3=1 



Qi(z,n) = 



qi 2 (z,n) 




Q 2 (z,n) = 





q 2 i (z,n) 



U 2 (z,n) 



J 

Pn+l 





satisfying the equations 

Uj(n) + A(z)Qj(z,n) -Qj(z,n + l)A(z) = 0, j = 1,2 (3.5) 

with 
respectively. 

Some useful properties of the above functions are 

Lemma 3.1. With the conditions given by (3.1): 

(1) For \z\ < 1 qij(z,n) — > as n — > oo, 6oi/i are /mite /or almost every 
z G T; 

(2) For almost every z G T (721 (£j n) — > as n — > 00; 

(3) For |z| < 1, qij(z,n) G ^2 and belong to H 2 as functions of z for all 
n > 1. 

Proof. For |z| < 1, it was noted by Benzaid and Lutz that (1) and (3) follow 
since qij is the convolution of an I 2 sequence with an l l sequence. For \Q\ = 1 
the results follow from the celebrated Carleson's theorem. □ 

For \z\ < 1 define Y(z,n) from 

*n(*) = (/ + Qi(z, n))(I + Q 2 (z, n))Y(z, n). (3.6) 

Then because of (3.2) and (/ + Qj)" 1 = I — Qj, we have 

Y(z, n + 1) =(I - Q 2 (z, n + 1))(J - Q^z, n + 1))(A + U(n))(I + Q^z, n)) 
(I + Q 2 (z,n))Y(z,n). 

After multiplying out the brackets the "first order terms" vanish (that is the 
key idea of Benzaid and Lutz) , and we come to 



Y(z, n + 1) = (A(z) + V(z, n))Y(z, n) 
with 



V(z,n) 



vn(z,n) v 12 (z,n) 
v 21 (z,n) v 22 (z,n) 



vu(z,n) = P n+1 q 21 (z,n)(l + p(z,n + 1)), 

v 12 (z, n) = f3 n+1 (p(z, n + 1) + p(z, n) + p(z, n)p(z, n + 1)), 

v 21 (z,n) = -/3 n+1 q 21 (z,n)q 21 (z,n+ 1), 

v 22 (z,n) = -/3 n+1 q 2 i(z,n + 1)(1 + p(z,n)), 



(3.7) 



where p(z,n) = qi 2 {z,n)q 2 \(z,n). By Lemma 3.1 V(z,n) G l\ for each 
\z\ < 1. 

Let T r = {\z\ = r} be a circle of a radius < r < 1. For z G T r we define 
W(*,n) = A- 1 F(z,n) G £\ then 



Y(z, n + 1) = A(z)(J + W(«, n))y(z, n). 



Consequently, for z £ T r 

\\Y(z,n + l)\\ < \\A(z)\\ \\I + W(z,n)\\\\Y(z,n)\\ 

< \\I + W(z,n)\\\\Y(z,n)\\ < (1 + || W(z, n)\\)\\Y{z, n) 

oo 

< ||Y(z,0)||exp^ \\W(z, n)|| < Cl , 

n=0 



(3.8) 



where 



ci= sup ||F(z,0)||exp]T||W(z,n)||. 



It now follows by induction that 

- z n-m o 



Y(z,n) - 



1 



y(z, m) 



< cl £pn^)ll- 

k=m 



Let us write the latter inequality for the vector function Y(z, n) 
coordinatewise 



(3.9) 

yi(z,n) 
V2(z,n) 



\yi(z,n) - z n m yi(z,m)\ + \y 2 (z,n) -y 2 (z,m)\ < c x sup ^ ||VF(z,/c)|| 

\ z \= r k= m 

(3.10) 

for \z\ = r, and hence by the Maximum Modulus Theorem for \z\ < r. So 
lim yi(z, n) = 0, lim j/2(-z, i) = u(z\ a, /3), (3-11) 

n— >oo n— >oo 

u is an analytic function, uniformly on compact subsets of the unit disk. 
From the relation between Y(n) and 4> n (3.6) 



1 qi 2 (z,n) 
q 2 i(z,n) l + p(z,n) 



yi(z,n) 
y2{z,n) 



(3.12) 



and Lemma 3.1 we see that the above implies that {<& n (z)} converge uni- 
formly on compact subsets of \z\ < 1 to 





u(z;a,P) 



(3.13) 



Since by (2.1) the leading coefficients K n converge, the same conclusion 
holds for 3> n . 

Now, let \z\ > 1. By using (2.6) we see that ^ n converge uniformly on 
compact subsets of the unit disk to 



Thus we have proved, 





v(z) 



v{z) = u{z 1 ;(3,a). 



(3.14) 
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Theorem 3.2. Suppose a n and f3 n satisfy (3.1). Then & n (z) converges 
uniformly on compact subsets of \z\ < 1 to (3.13), and ^f n (z) converges 
uniformly on compact subsets of \z\ > 1 to (3.14). 



4. Tauberian results 



We now consider what happens for \z\ = 1, and also the boundary values 
of the functions considered in the previous section. Although {z 11 } is no 
longer in £ l , the above argument goes through as long as we require further 
assumptions on the coefficients. Similarly to (3.4) we define 



qn(z,n) = ^a k z k n , q 21 ( 



z,n) 



k=l 



k=n+l 



n—k 



(4.1) 



Clearly, Lemma 3.1 holds for q~ij as well if \z\ < 1 is replaced by \z\ > 1. 
Next, put 



M = < z = e : 



k=i 



+ 



E 

k=l 



a k z 



< oo 



(4.2) 



and 



E = {z = e w : J2 \Pk+i<i2i(z, k)\ < oo}, 



fc=o 

oo 



(4.3) 



E = {z = e ie : \ak+i<hi(z, k)\ < oo}. 
fc=0 



We know from the theory of Fourier series that under assumption (3.1) M 
is the set of full Lebesgue measure, and for z 6 M 

oo oo 

qn(z, n) - q2i(z, n) = z n ^ P k z~ k , q 12 (z, n) - q 2 i(z, n) = z~ n ^ a k z k . 

k=l k=l 

Theorem 4.1. Suppose (3.1) holds, and z G E PI M, then for u n (2.2) 

lim u n (e id ) =u*(e id ) 

exists. Likewise if (3.1) holds, and z G E D M, then for v n (2.3) 

lim v n {e ie ) =v*(e td ) 

n->oo 

exists. 

Proof. For z € M we have 

sup \q\ 2 (z,n)\ = C(z) < oo, 



so for the matrix entries Vij (3.7) 

|t>n (2, n) | < (1 + C\q 21 (z,n)\)\l3 n+1 q 21 (z,n)\, 

\v 12 (z,n)\ < C\(3 n+ i\(\q 2 i(z,n) \ + \q 2 i(z,n+ 1)| + C\q 21 (z,n)q 21 (z,n + 1)|), 

\v 21 (z,n)\ < \f3 n+1 q 21 (z,n)q 21 (z,n + 

|t> 22 (z,n)| < (1 + C\q 21 (z,n)\)\(3 n+1 q 21 (z,n + 1)|. 

(4.4) 

Also 

\z[3 n+1 q 21 (z,n + 1)| < |/S n +i92i(^,n)| + l^+i^i^n + 1) - g 2 i(^,n))| 
= |/S n +i92i(n)| + |/3„+ia n+ i|. 
Thus if z € -E 1 fl M, we find 

00 

^||F(z,n)|| <oo, 

n=0 

and as in (3.8) above 

sup||F(z,n)|| = Ci(z) < 00, z€ fin M. 

As in (3.10) we see that 

n 

\y 2 (z,n)-y 2 (z,m)\<C 2 J2\\V(z,k)\\, 

k=m 

which shows that linin^oo y 2 (z, n) = u(z;a,/3) exists and is finite for z £ 
EDM. From the relation between $ n and Y 2 (3.12) and (2.2) we derive 

K~ 2 u n (z) = (1 + p(z, n))y 2 (z, n) + q 21 (z, n)yi(z, n). 

Since yi(z, n) = 0(1), n ->■ 00 (cf. (3.10)), and 921(2, ra) -> for z G £ n M, 
we find 

00 

u*(z) = lim u n (z) = Y\(1 - aj^y 1 u(z;a,f3). 
i=i 

An analogous argument shows that if z € E n M then 

00 

u*(z) = lim f n (z) = TT(1 — aj/3,-) -1 ^(z" 1 ; /3, a) 

exists and is finite. □ 
Next, set 

00 

E{6)= sup X> fc+ ig 2 i(re* fc)|, 



°<^ 1 'fc=0 



(4.5) 



£(0) = sup £|a fc+ i«ai(re* fc)l> 



l<r<oo fc=0 
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and 

N(9)= sup \q 12 (re ie ,n)\, N(0) = sup \q 12 (re ie ,n)\. (4.6) 

0<r<l,n l<r<oo,n 

We now examine what happens when z = re td approaches the unit circle for 
certain values of the argument 8. 

Theorem 4.2. Suppose (3.1) holds, and z = e ld is such that E(8) < oo and 
N(8) < oo, then for u given by (3.11) the limit 

lim u{re ie ) = u*(e ie ) 

r— s-1— 

exists. Likewise if z = e ld is such that E{8) < oo and N(8) < oo, then the 
limit 

lim v(re ld ) = v*(e td ) 

exists. 

Proof. Suppose that the first assumption holds. This implies that inequali- 
ties (4.4) hold with C replaced by N(8). Since E(8) < oo we find that 



sup \^ ||V(re i6, ,n)|| < oo. 
Thus sup 0<J .< 1 \\Y(re ie ,n)\\ < cN{8) and 

n 

sup \y 2 (re te ,n)-y 2 (re t8 ,m)\ <cN(8) sup V \\V{re w , i)\\. 

0<r<l 0<^<1~; 

Hence y 2 (re te , n) converges uniformly for < r < 1 to u. So u n (z) converges 
u*(z) uniformly for < r < 1. Since u n are continuous for < r < 1 (they 
are polynomials) we find that lim r _>i_o u(re 10 ) = u*(e %e ). An analogous 
argument proves the second statement. □ 

Finally, we conclude that if z G M, and both E(8) and N(8) are finite, 
then 

lim lim u n (re tS ) = lim lim u n {re lS ). 

r->l-0n->oo n— >oo r— ¥l— 

Similarly, if z G M, and both E(8) and N(8) are finite, then 
lim lim v n {re 1 ®) = lim lim v n (re i ). 

) — ^1+0 n— too n— <roo r— >l+0 

In the setting of polynomials orthogonal with respect to a positive measure 
fi on the unit circle ([3 n = a n , \a n \ < 1), the tauberian problem in the Szego 
class was studied first in [8, Chapter 5]. In this case M = M and 

oo oo 

E = E = {z = e ie : ^ \a k +i Yl a j zj ~ k \ < °°}- ( 4 - 7 ) 

fc=0 j=k+l 

Condition (4.7) appeared in [14], where it was shown to imply a bound 
for orthogonal polynomials | log |$ n (z)|| = O(l) (see also [12]). (4.7) is a key 
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ingredient in [11], where measures on the unit circle with slowly decaying 
parameters are studied and the Benzaid-Lutz method is applied. It is shown 
there that the uniform convergence in (4.2) on an arc along with (4.7) implies 
the uniform convergence of the reversed polynomials on the same arc (p, is 
said to admit a uniform asymptotic representation). Damanik [5] proved 
that the singular component fx is supported on the complement of E. 

We complete with two examples (the first one is borrowed from [17, Chap- 
ter V.4]). 

Example 4.3. It is known that the Fourier series 

00 ick log k 

^) = E^7^ e ^' 0<a<1 ' c>0 

fc=l 

converges uniformly on T, and g a £ Lip(a). Hence (see [17, formula 
(13.26)]) 

q 2 i(e %e ,n) = 0(n~ a log n) 
uniformly on T, and the measure fi with parameters a n = n - 3 / 4 - £ e * cnlo s n ) 
e > 0, admits the uniform asymptotic representation on the whole circle. 
The same result holds for fi with a n = n~ 3 / 4-e e ma , < a < 1. 

Example 4.4. Let 

1 m 

3=1 

for large enough n > n . We have 

oo m oo ^ 

-e-V(^,n)= «^ = £^ E y^+V 

k=n+l j=l k=n+l 

so the series converges uniformly inside T\{£i, . . . , Cm}, Ck = e~ tXk with 
the bound 0(l/n 7 ). Hence for 7 > 1/2 the measure with such parameters 
admits the uniform asymptotic representation on T\{Ci, . . . , Cm}- 
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